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ABSTRACT 

The field of extrasolar planets has rapidly expanded to include the detection 
of planets with masses smaller than that of Uranus. Many of these are expected 
to have little or no hydrogen and helium gas and we might find Earth analogs 
among them. In this paper we describe our detailed interior models for a rich 
variety of such massive terrestrial and ocean planets in the l-to-10 M e range 
(super-Earths). The grid presented here allows the characterization of the bulk 
composition of super-Earths detected in transit and with a measured mass. We 
show that, on average, planet radius measurements to better than 5%, combined 
with mass measurements to better than 10% would permit us to distinguish 
between an icy or rocky composition. This is due to the fact that there is a 
maximum radius a rocky terrestrial planet may achieve for a given mass. Any 
value of the radius above this maximum terrestrial radius implies that the planet 
contains a large (> 10%) amount of water (ocean planet). 

Subject headings: planetary systems — planets and satellites — Earth 
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Introduction 



A new family of exo-planets is now being explored thanks to improvements in detection 
methods that allow the discovery of very low mass objects (1-10 M @ hereinafter super- 
Earths). The first three super-Earth s discovered are: GJ876d - a 7.5 M® planet at 0.02 
AU from its star (iRiyeraelalJ 5)05]) , OGLE-2005-BLG-39 0Lb - a 5 M m pla net at 5 AU 



(feeaulieu et aDbood ) and HD69830b - a 10.8M e at 0.08 AU Lovis et a 



2006). In the next 



few years, with missions like COROT, that launched in Decemb er 2006 (IBorde et al.l (120031 )) 
and Kepler, scheduled to launch in 2008 (IBorucki et al.ll2003l ). many more - and smaller, 
super-Earths are expected to be found. These missions detect planets by observing their 
transits, hence allow for a measurement of the planet's radius, R. The Doppler techniques 
have improved significantly and allow detection and planet mass, M, determination for 
super-Earths as well, as demonstrated by the discovery of GJ876d and HD69830b. 

Our work is motivated by the need for mass-radius relations for super- Earths of varied 
compositions. Such relations should allow discerning the bulk compositi ons of planets de- 
tected in transit - with derived R and M. Ze ro-temperature sphere models (jZapolsky fe Salpeter 
19691 ; IStevensonl 11982c iFortney et al.l 120061 ) allow approximate mass-radius ca l culati ons for 
specific simple materials (pure Fe, pure olivine, pure water). IValencia et al.l (120061 ) intro- 
duced detailed interior mo dels for super-Earths. These models were extended to sup er-Earths 



with l a rge water content ([Valencia et al.l 120071 ) - the ocean planets prop osed by 



Kuchner 



(2003); iLeger et al.l (12004 ). A cold ocean planet model was introduced by lEhrenreich et al 
d2006f ) for the low-temperature planet OGLE-2005-BLG-390Lb. 



Distinguishing between ocean planets and dry rocky super-Earth planets is not trivial, 
because at a given mass the sa me planet radius coul d correspond to different mixtures of iron 
core, mantle minerals, and ices ([Valencia et al.ll2007l ). Observations with certain t olerances of 



the de rived planetary R and M are required, as shown for two particular cases by lSelsis et al. 



( 120071 ). On the other hand, there is very strong motivation to be able to distinguish between 



ocean planets and dry planets, because of all the implications to planet format ion theory 



2004; 


Rafikov 


2006: 


Kennedv et al. 


2006; 


Ravmond et al. 


2006 



anticipate the formation of 
super-Earths with a rich variety of bulk compositions, but differ in the details. 

As the field of planets with masses in the range of 1-10M® emerges, it will be useful to 
have a consistent nomenclature to describe such planets. We propose to call a super-Earth, a 
planet with large amounts of solid material in the l-lOM® range. The lower bound is obvious 
for historical reasons, the upper bound is somewhat arbitrary. It is informed by the physical 
argument that a t ~ 10M ffi and above, the planet can retain H and He during formation 
( llda fe Linll2004l ). Super-Earths, includes ocean planets and massive terrestrial planets as a 
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subclass in this mass range. Planets with H2O on the order of > 10% in either liquid and/or 
solid phase are termed ocean planets. Massive terrestrial planets will then be the ones with 
little water on the surface or in the mineral assemblages and hence, be the most similar to 
Earth. The existence of a large core of iron makes a difference, therefore within super-Earths 
there is a subset of planets with large cores that we have termed 'super-Mercuries', and which 
can be terrestrial or have a large ocean. The name 'ice giants' - Uranus and Neptune - has 
been reserved for planets more massive than 10M ffi with some amount of gas (~ 10%). In 
addition, it is possible, though unlikely, that planets with M < 10M ffi could have a lot of 
gas (> 10%), but we do not model them in this paper. At 10M ffi , a complication arises due 
to the possibility of H and He retention by the planet during formation. We do not consider 
these planets in this paper as a separate class since they can be treated as massive ocean 
planets with high H 2 content and a relatively small layer of gas above it that would add a 
few kilometers to the radius. The regime of M < 1M®, that includes Mars, Venus, Mercury 
and the icy-satellites of the Solar System has been studied previously with models describing 
these planets and therefore has been left out of this paper. 

In this study we explain how we can differentiate between a rocky and an ocean super- 
Earth with measurements of mass (M) and radius (R) at the current precision levels. In the 
first section we briefly explain the model used to obtain the internal structure of a super- 
Earth given its composition (%) and M. This model has been used to determine the total 



k parameters 



radius and interior distribution of bu 
of terrestrial and ocean exo-planets (Valencia et al. 



mass, gravi ty, pressure and temperature) 



20061 . 120071 ) and here we present another 



useful application of this model: how to obtain constraints on \ with values of M and R. 



Model 



We model the internal structure of a low mass exo-planet by dividing its interior in con- 
centric shells where composition is considered homogenous. Within each layer the numerical 
model solves the standard differential equations for density, pressure, mass and gravity struc- 
ture under hydrostatic equilibrium. The equ ation of state (EO S) chosen to relate density, 
pressure and temperature is the Vinet EOS (IVinet et al.l Il989f ) that e xtrapolates better to 



high pressures th an the third-order Birch-Murnaghan (BM3) EOS (IHama fc Suitol 11996 
Cohen et al.ll2000l ). The regions are: 



an H 2 layer divided into a water ocean over an icy shell - composed of high-pressure 
ices VII and X 



the mantle divided into an upper mantle - composed of olivine ([Mgi-^Fe^SiO^ 
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wadsleyite ([Mgi-^Fe^SiOzi)^ and ringwoodite ([Mgi^Fe^SiCU) 7 - and a lower 
mantle - composed of perovskite ([Mg 1 _ a .,Fe a! ]Si03), ferromagnesiowusite ([Mgi^^e^jO) 
and post-perovskite ([Mgi^^e^jSiOs) 

• and the core (mainly Fe), that can be divided into a liquid outer layer and inner solid 
layer if the temperature profile of the planet crosses the melting curve of iron. 



The thickness of the H2O, mantle and core regions are dependent on the ice mass fraction 
(IMF) and core mass fraction (CMF) that are assumed for a planet. The Earth has a 
negligible wat er fraction (IMF=0.02-0.1%) and a core that is 1/3 of its mass (CMF=32.59%, 
Staceyl (11992 )). Phase transitions are modeled in each region to obtain the different layers 
composing the planet. The temperature profile is modeled after the Earth's; the temperature 
gradient is adiabatic within each region and conductive within the boundary layers developed 
at the top and bottom of the mantle. The boundary layers in the H2O and core regions are 
neglected owing to their very low viscosity that prevents thick boundary layers to develop. 
The surface temperature is taken to be 500 K after the d etection of GJ876d with an estimated 
surface temperature of 430 — 650K (IRivera et al.l 120051 ) . Despite this high value, water may 



be present provided the surface pressure is large enough f ~ 1000 atm for T sur f = 500). 
For more details please refer to IValencia et al.l (120061 . 120071 ) where the model is explained 



thoroughly. Currently the models are not evolved in time; they correspond to a planet that 
has undergone differentiation and achieved equilibrium like our Earth. 

This numerical model allows us to calculate the total radius of a planet given its mass and 
composition. We do not know the exact composition of the planets that will be discovered 
in the near future but we can make reasonable assumptions on their bulk composition. The 
dominant minerals in the mantle are most likely the same as for Earth (silicates with a small 
incorporation of iron -~10%). The core is mainly composed of Fe. The Earth's core is known 
to have some Ni (mostly in the inner core) and a light alloy in the outer core of ~10%. The 
nature of the light alloy and exact amount are not known for Earth. Instead of looking at all 
possible candidates (O, H, S, Si, C) we model the core as pure Fe in this study and anticipate 
that the presence of a light allow will increase the radius by a few hundred kilometers. The 
water region is composed of ice VII and X. For low temperature planets, there might also 
be an extra icy layer of Ice I above the water ocean. This layer will be small due to the 
small stability field of Ice I and hence we consider the results presented in this study to 
be adequate even for low temperature icy planets (such as Ganymede, Europa and Calisto 
massive analogs). For planets that are very close to its parent star, their surface temperature 
might be on the order of 500 — 1500K. Ice VII is a high pressure form of H2O that could be 
the source of a liquid water la yer, provided the surfac e temperature is high enough according 
to the phase diagram of H 2 (jWagner &: Prussll2002l ). The thermal structure has little effect 
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on the radius of a planet (I Valencia et al.ll2006l ) but imposes conditions on the compositional 
phases, particularily of H 2 at the surface. If the surface temperature is higher than the 
critical point (T crit = 647K and P cr it = 218 atm) the structure of an ocean planet will vary 
gradually from a vapor phase to an ice VII phase at depth. 



3. Ternary Diagram 

We have computed a grid of interior models appropriate for ocean and terrestrial planets 
with the entire range of CMF and IMF to determine how M, R and x are related. Since 
the percentage of mass in the mantle, ices and core has to add up to 1, we use a ternary 
diagram to show the results. Ternary diagrams are commonly used in Earth sciences. Data 
of a three-component system is plotted in a triangle whose sides depict three axes. Each 
vertex means 100% of a particular component and data plotted on the opposite side means 
0%. Lines parallel to a particular side will show various degrees of a component whose 
maximum value is shown at the corresponding opposite vertex. We use a fourth axis (color 
code or label) to show the value for the radius for each composition. A 3-dimensional volume 
describes the complete family of Mass-Radius relations for all compositional combinations; 
a ternary diagram is a cross-section at a given mass. Figure 1 shows some examples for a 
5M e planet (like OGLE-2005-BLG-390Lb) with three different compositions (A) CMF=10% 
and IMF=50% with a radius of 12200km, (B) CMF=20% and IMF=10% with a radius of 
10750km , (C) CMF=30% and no water with a radius of 9800km and (D) CMF=40% and 
IMF=30% with a radius of 10950km. 

A minimum and maximum radius exists for a planet with mass M. The smallest a 
planet can be is when it is composed of pure iron or heavy iron alloys (Fe+Ni) and is shown 
in the right vertex of the ternary diagram. The maximum size for a semi-solid planet (not 
gaseous) is obtained if it is composed purely of ices and water and plots in the left corner. 
For a lM e planet these extreme values are ~4900 km and ~9300 km. Simple formation 
constraints suggest that the minimum and maximum are not quite at the vertices (more 
realistic conditions see section 4.1). Terrestrial planets (no substantial water content) of 
5M® are shown in Fig. 2 on the side connecting mantle and core only. 
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4. Results 

We have calculated the ternary diagrams for planets with 1M®, 2.5M e , 5M ffi , 7.5M e 
and 10M® showing the maximum and minimum radius for each planet and curves of con- 
stant radius. Figure 2 shows the different ternary diagrams. The same value for radius 
can be obtained through different combinations of core, mantle and H 2 contents showing 
degeneracy in composition. 



4.1. Initial Conditions and Constraints on Bulk Compositions 

Our interior model of a super-Earth planet could not be a random mixture of layers of 
different mineral, alloy, or compositional shells. This is so for at least two general reasons: 
differentiation and constraints in the proto-planetary disk abundances. Therefore we make 
two basic assumptions in this work. The first is that all our model planets have undergone 
differentiation. The second is that the relative elemental abundances are the same as those in 
the Solar System and the solar neighborhood. The first assumption is not restrictive since all 
the terrestrial planets and large satellites in the Solar System are known to be differentiated. 
The second assumption should also be adequate for the super-Earths that we expect to be 
discovered in the near future (close-by and in the pre-selected targets of the Kepler mission). 

Thus, we consider three possible initial states for our models: 

(1) formation and post-differentiation anywhere in the proto-planetary disk - the ratios 
of Si/Fe, O/Fe, etc. depend on T con d and the amount of Fe incorporated in the mantle 
depends on the O locally available; 

(2) formation as in (1), followed by a collisional stripping - then four initial states are 
possible: (a) pure iron alloy core; (b) core and mantle; (c) pure mantle; and (d) mantle and 
water ocean; and, 

(3) formation as in (1) or (2), followed by water delivery - then two initial states are 
possible: (a) core, mantle, and water ocean, and (b) mantle and water ocean. 



To explore the first scenario we consider that during formation, the building blocks of 
planets - volatiles, silicates and metals - condense out at different temperatures T con d. As 
the nebula cools down the most refractory elements will condense out first. The first of the 
building blocks of planets to condense out are silicates at temperatures bet ween 1750-1060K 
fcrot et al.ll2005h . followed by the metals (Fe, Ni) between 1450-1050K fcrot et alJl2005h 
and finally H2O and other ices, depending on the pressure. The relative amounts of the 
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different major elements in solid planets are H at 74%, O at 1.07%, Fe at 0.1%, Si at 0.065% 
and Mg at 0.058%. We consider Si and Mg to be practically equally abundant. During the 
condensation sequence (for pressure < 10 _4 bars), Si will condense before Fe. If Fe remains 
immiscible it will form the core during condensation and planet segregation. This is the end 
state scenario that yields the largest core possible for a planet. That is, a fixed ratio of Si/Fe 
of ~ 0.6. The mantle is effectively MgO+Si02 so that Si/Fe ratio can be used as a proxy 
for MMF/CMF. This line is shown in Fig. 1 as the lowest possible values for MMF/CMF. 
Any value above this line is plausible during formation because Fe can be incorporated in 
the mantle to form (Mgi_ x , Fe)0+SiC>2 (with x = 0.1 for Earth) depending on the redox 
conditions. 

To explore the possibility of having a planet with a large core and H 2 layers and a 
small mantle we consider a scenario with a pure Fe embrio that acquires water from comet 
delivery, which is scenario (3) stated above. Comets have a dust/gas ratio of 1-2. As 
the comet delivers water (gas, volatile) it will at least deliver Si+Mg+Fe+O in the same 
proportions. We consider that the cometary dust is made up of Si, Mg, Fe and O and that 
the proportions of these elements are the same as in the solar nebula (or CI chondrites). 
Then the ratio of Si/H 2 that needs to be maintained is at least 0.23 by mass. Therefore, 
for every gram of water delivered to the Fe pure embrio, 0.23g of Si (that makes up the 
mantle) will be delivered also. This is the end state case that yields the lowest amount of 
MMF possible for a given IMF. This line shows the minimum amount of mantle formed as 
planetesimals get accreted and deliver water. Accompanying this amount of Si, is the amount 
of Fe that, if in solar abundances can also add mass to the core. Therefore the unshaded 
region in Fig. 1 shows the plausible compositions that reflect the solar nebula composition 
and cometary delivery constraints. 

To explore scenario (2) mentioned above, we contemplate late-stage processes that might 
influence the final accretionary state of a planet starting from a 'normal' embrio. Any process 
that induces the escape of preferentially light elements (solar wind, gravitational escape) will 
deplete the planet from volatiles, H2O and perhaps silicates forming terrestrial planets (core, 
mantle+core planets). Conversely, any process that targets heavy elements to be excluded 
from accretio n (like formation of the Moon after the Mars-sized impact onto a proto-Earth 



(ICanupll2004j )) would deplete the planet from Fe and silicates forming a light planet (mantle, 
mantle+ water ocean planets). It is difficult to imagine a process that would preferentially 
retain the heavy and light elements simultaneously and blow-off or prevent from accreting 
the medium density components (silicates). 

Based on these arguments we expect that planets with large CMF and IMF and small 
or non-existent mantles are very unlikely. Thus, the shaded region in Fig 1 is unlikely to be 
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populated. 



4.2. Maximum Radius for Terrestrial Planets 

For a given planet, a threshold in radius exists such that larger values imply that the 
planet i s ocean-like - has a su bstantial amount of water and ices and can not be completely 



rocky ( (I Valencia et al.l 120071 ). This threshold radius can be identified as the largest iso- 
radius curve that intersects the terrestrial side on Fig 2. All the iso-curves to the left (bluer 
region) of this curve show IMF>0. This threshold in radius is 6600 km, 8600km, 10400 km, 
11600km and 12200km for lM e , 2.5M e , 5M ffi , 7.5M e and 10M ffi planets respectively Thus, 
for a given planetary mass, a measurement in radius larger than its corresponding terrestrial 
threshold radius would imply that the planet either formed beyond the snow line acquiring 
large amounts of H2O, or had large amounts of volatiles delivered to it by wet planetesimals. 

The identification of the composition via this maximum terrestrial radius will be easier 
for smaller planets since their iso-curves are more significantly separated in the ternary dia- 
gram. 



4.3. Iso-radius curves 



Figure 2 shows the progression in radius from a 10M e to a 1M® planet and how the 
iso-radius curves move closer to each other as the mass increases. The lines of constant 
radius are less sensitive to the amount of mantle mass, illustrated in how parallel they lay 
with respect to increasing values of mantle percentage (iso-curves almost vertical). 

• This is particularity true for a 5M e planet: a large range of mantle mass (between 
— 85%) may yield the same radius of ~ 11300km, provided the correct amount of 
water is present (from 50% to 15% respectively). 

• For a lM e planet, the iso-radius curves are tilted to be slightly more sensitive to 
mantle mass, less sensitive to the amount of core and more sensitive to water content 
(iso-cuves are more parallel to the lines of constant water content). At this pressure 
regime (smallest planet shown), the amount of low density (water) and high density 
(Fe) material can proportionally affect the total radius. 

• For the case of a 10M e (high internal pressures) the iso-curves show that for very 
water rich planets ( > 70% H 2 0), the radius is dependent on the amount of core-mass 
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fraction (blue lines moving towards being parallel to the terrestrial line) and become 
less independent as the composition becomes iron rich (> 70% Core). For the large 
core fractions, the amount of water is important because the contrast in density is the 
largest. The overall trend is that the iso-radius curves are less sensitive to the amount 
of mass in silicate form (mantle) and the result for radius depends more in the amount 
of light elements (water and ices) and heavy elements (core-Fe) . 



There h ave been useful mass-radius relationships obtained with zero-te mperature sphere 
models (jZapolsky fe Salpeterl 1 19691 : IStevensonl Il982l ; iFortney et al.l 120061 ) that hint to the 
bulk composition of a planet. These models find the radii for different pure compositions 
such as H2O, silicates or Fe. These radii (as a function of mass) correspond to each vertex 
of our ternary diagram. They also mix pairs of the above compositions, essentially following 
the sides of the ternary diagram. It can be seen from Fig. 2 that the iso-radius curves are 
very sensitive to the mixture and how they intersect with the ternary sides depends on it, 
especially in a mantle-H^O mixture. Of course, the zero-temperature sphere models have 
no layering, or phase transitions, or convective mixing in the planet's interior. They also do 
not access the parameter space inside the ternary diagram. Therefore, a zero-temperature 
model is a good first approximation. However, detailed calculations are needed to distinguish 
between realistic compositions at the precision C0R0T and Kepler will allow. 



5. Mass-Radius Relationship 



We had previously stated that the radius of a super-Earth scales in a power law rela- 



tionship with mass: R = aR @ yjj^j where a = 1 and (3 ~ 0.27. We have improved this 
result by choosing an equation of state better suited for high-pressure extrapolations (Vinet 
EOS) and including the rec ently discovered new si licate phase - post-perovskite (ppv) in the 
Earth's lower-most mantle (IMurakami et al.ll2004l ). Even though the stability field of ppv is 
only of ~ lOGPa on Earth, it constitutes most of the mantle of super-Earths. It also can 



accomodate more Fe in its structure than perovskite (IMao et al.ll2004l ) making it more dense. 
With these improvements the exponent of the power law becomes smaller: f3 = 0.262. It is 
different from 1/3 - the constant density scaling value - due mostly to pressure effects. The 
scaling is independent of temperature effects (such as surface temperature, constant versus 
temperature-dependent viscosity cases), and CMF. It also differs little with core composition. 



Sotin et. al (jSotin et al.l 120071 ) has shown that planets that have 50% water content 

/ \ 0.274 

scale as R = 1.25 x R m ( 4r- I for a family of planets that have the same Fe/Si ratio. 



Al: 



Here, we generalize the power law relation for low-mass planets to include up to 50% H 2 
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assuming a fixed ratio between the silicate and metal portion of a planet. Meaning, that the 
mantle and core mass fractions are the same so that the Fe/Si ratio of the family of planets 
with different IMF varies little. Figure 3 shows the radius as a function of mass in terms of 
Earth-masses for rocky and ocean planets. The generalization of the power law becomes 

i2p = (1 + 0.56 x IMF) i2e(-/ 

where IMF denotes the amount of water in percentage. Table 1 shows the values for the 
coefficient and the exponent in the power law corresponding to the family of planets with 
up to 50% H 2 0. The coefficient a shows how much larger Earth would be if it had different 
amounts of water. Earth would have an increase of 26% in radius if it had 50% water by 
mass. Our result for a 50% H 2 0-planet differ from the exponent obtained by Sotin et. al. 
possibly because we in clude the newly discovered post-perovskite phase in the lower mantle 



([Murakami et al.ll2004j ) and our preference for the Vinet EOS over the BM3 EOS chosen by 



Sotin et. al. We are presently working towards an understanding of this difference. 

Figure 3 shows the mass-radius relationship for planets with a fixed mantle-to-core ratio 
of 2:1 as a proxy for Fe/Si. It also shows the relationship between the planetary mass and 
both the minimum and maximum radii. The minimum radius is obtained for a rocky planet 
with the largest possible core constrained from solar nebula composition. The maximum 
radius occurs when the planet is made of silicates and water constrained by the volatiles/Si 
ratio from chondritic material (refer to section 4.1). 



6. Observational Characterization of Super-Earth Planets 

The main result from our model calculations presented here is that there is a threshold 
of precision in R and M that allows bulk estimates of the composition \. These are ~5% 
in planetary radius and ~ 10% in planetary mass (this is an average statement over the 
mass range 1 — 10M ffi and our complete grid of models). Instrumentation currently under 
development will be able to provide such observational tolerances in the coming 5 years - in 
particular, the NASA Kepler mission in synergy with the HARPS-NEF spectrograph. 



The Kepler mission (IBorucki et al.ll2003l ) will monitor a single field in Cygnus/Lyra and 



discover many hundreds of transiting planets orbiting preselected main-sequence stars (F to 
M-type) of 11th to 15th magnitude in V. It is capable of detecting more than 50 Earth-sized 
(li?e) planets in all orbits, with about a dozen n orbits approaching 1-year period; as well 
as several hundred super-Earths (~ 1.3i?®) in all orbits up to 1 year. The Kepler mission 
will achieve 20 ppm photometry or better on 12th magnitude stars (its "sweet spot" for 
detection). 
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The HARPS-NEF spectrograph is a collaborative Harvard-Geneva Observatory project 
to build and operate a high-precision radial velocity instrument in the Northern h emisphere 
for synergy with Kepler. Its design is based on that of HARPS (jPepe et al.ll2005l ); it will be 
capable of achieving better than lm/s in 1 hour on a 12th magnitude star of F- to M-type. 
That would allow planet mass determination of 10% or better for super-Earths (2-5 ®) in 
orbits closer than 0.10 AU. For more massive super-Earths, 10% in mass is achieved for 
orbits up to 0.5 AU, all for planets already detected with transits by Kepler. 

Using Kepler and HARPS-NEF together should allow tolerances in radius (< 5%) 
and mass (< 10%) for dozens of super-Earths, limited only by available observing time for 
HARPS-NEF. This would produce a large enough statistical sample to distinguish the bulk 
properties of ocean planets and dry rocky planets and their rate of occurance in orbits closer 
than ~1AU. Examples of the ability to distinguish these properties are given in Fig.4 and 5. 

Our estimate of 4% uncertainty in planet radius from Kepler incorporates the photo- 
metric uncertainty for 10 high-impact transits on a 12th magnitude star and the uncertainty 
in the stellar parameters (jCody fc Sasselovll2002l ). and stellar radius in particular. This es- 
timate makes use of the fact that Kepler will discover a large number of planets and only 
a sub-sample of the most suitable ones will be followed-up with HARPS-NEF and will have 
Kepler-derived parallaxes. A g ood general discuss ion of observational uncertainties for both 
CoRoT and Kepler is given by ISelsis et al.l (120071 ). 



The apparent radius of a transiting planet depends on its atmosphere; iBurrows et al. 



( 120031 ) have estimated a several percent increase in the radius of a hot Jupiter compared to 
a reference radius at 1 bar pressure. We have neglected this 'transit radius' effect here, but 
note that it could be important for super-Earths with significant atmospheres. 



7. Conclusions 

We have present here the relationship between planetary mass, radius and bulk com- 
position for planets in the 1M® to 10M© range. The main bulk composition materials can 
be described well in three groups: core (Fe), mantle (silicates) and water (H2O). Therefore 
the planet radius dependence on bulk composition for given mass is best illustrated with the 
help of a ternary diagram. The degeneracy in composition is illustrated in curves of constant 
radius. These iso-radius curves are less sensitive to mantle composition, showing that the 
degeneracy is mostly achieved by a trade-off between core and H 2 content. 

Our detailed models indicate that a maximum terrestrial radius exists for a given plan- 
etary mass - it separates dry rocky compositions from ones containing 10% or more water 
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by mass. A measured planet radius that exceeds that critical terrestrial radius would imply 
that the planet is an ocean planet. For the high-mass range, the possibility of a massive 
hydrogen-helium envelope has to be evaluated separately. 

By considering realistic initial states for our interior models we can constrain the pa- 
rameter space of plausible evolved models, as shown in Figs.l, 4, and 5. These constraints 
lead to minimum and maximum radius as a function of planet mass, and show that planets 
with small mantles and large core and H 2 contents are unlikely to form. 

We provide a mass-radius relationship to include the effects of various amounts of water 
and note that the exponent (3 in R = ai? e (M/M®)' 3 decreases slightly as the amount of H 2 
increases (due to the high compression of H 2 0) and a increases almost linearly with H 2 
content (due to the low density of H 2 0). 
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Table I. Generalized mass-radius relationship - R — aR 9 (-jj-) 



H 2 


a 


P 





1 


0.262 


10 


1.076 


0.260 


20 


1.130 


0.257 


30 


1.178 


0.252 


40 


1.223 


0.248 


50 


1.265 


0.244 


Note. 




Val- 



ues in the power law 
relationship between 
mass and radius of a 
super-Earth with sim- 
ilar Fe/Si ratio to 
Earth's 
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Fig. 1. — Ternary Diagram. The possible compositions of a super-Earth and radius are 
illustrated in a ternary diagram that depicts three axis: amount of H 2 0, core (Fe) and 
mantle (silicates). Each point inside the triangle uniquely specifies a bulk composition and 
the corresponding radius - a few examples are labeled. The terrestrial side is the location 
where rocky planets would be shown. The minimum plausible radius is achieved with the 
largest core possible from solar nebula composition constraints and no water. The maximum 
plausible radius is shown for a planet with no core and maximum amount of water under the 
solar nebula composition constraints. The shaded region shows improbably compositions 
from solar nebula abundances and accretion processes. 
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Fig. 2. — Ternary Diagram relating M, R, and x- The iso-radius curves for possible bulk 
compositions are shown for a 10M®, 7.5M e , 5M e , 2.5M e and 1M®. As mass increases, the 
iso-radius curves shift to the right and rotate slightly. The progression from a 5000km (1M® 
Fe-) planet to a 16000km (10M e water-) planet is illustrated in the color bar. White regions 
on the ternary diagrams are areas that pose a numerical challenge on the model. 




Fig. 3. — Mass- Radius relationship for Ocean and Rocky Massive Planets. The solid black 
line is the power law relationship for completely terrestrial planets with l-10M e . The dashed 
lines above progressively represent the relationship for planets with 10%, 20%, 30%, 40% 
and 50% H2O. This family of planets have a fixed mantle to core proportion of 2:1. The 
minimum and maximum planetary radius relations with mass are shown as dotted lines. 
Venus, Earth, Uranus and Neptune are shown for reference. 
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Fig. 4. — Kepler and HARPS-NEF capabilities for a 2.5 M e planet. The shaded area around 
the terrestrial radius threshold illustrates the uncertainty from a 4% error in radius mea- 
surements and 5% (short-dashed line) and 10% (long-dashed line) error mass measurements. 
The range of possible interior structures (bulk compositions) is dominated by the radius 
uncertainty in this range of masses. The two dashed lines define interior structures excluded 
by the initial conditions - see Fig.l. 
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Fig. 5. — The shaded area around the terrestrial threshold radius illustrates a 4% error in 
radius measurement for a 7.5M e planet. The effect of uncertainty in planet mass measure- 
ment is illustrated in the bottom panel: the same radius value of 12000 km could correspond 
to very different compositions, if the planet mass is not known to better than 30%. In such 
a case it is not possible to determine whether the planet is an ocean planet or a dry rocky 
planet. The two dashed lines define interior structures excluded by the initial conditions - 
see Fig.l. 
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